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Pulsar electrodynamics has been built up by taking ingredients from two 
'f^/ '\ models, the vacuum-dipole model, which ignores the magnetosphere but includes 

C/^ ! the inductive electric field due to the obliquely rotating magnetic dipole, and the 

^ \ corotating-magnetosphere model, which neglects the vacuum inductive electric 

field and assumes a corotating magnetosphere. We argue that the inductive 
field can be neglected only if it is screened by a current, Jsc, which we calculate 



Ph' 



O 



r- 



c/2 \ for a rigidly rotating magnetosphere. Screening of the parallel component of 

the inductive field can be effective, but the perpendicular component cannot be 
^ \ screened in a pulsar magnetosphere. The incompletely screened inductive electric 

^ \ field has not been included in any model for a pulsar magnetosphere, and taking 

Q \ it into account has important implications. One effect is that it implies that 

the magnetosphere cannot be corotating, and we suggest that drift relative to 
corotation offers a natural explanation for the drifting of subpulses. A second 
O \ effect is that this screening of the parallel inductive electric field must break 

down in the outer magnetosphere, and this offers a natural explanation for the 
acceleration of the electrons that produce pulsed gamma-ray emission. 

X 

H \ Subject headings: stars: pulsars: general radiation mechanisms: non-thermal - 

magnetic fields 



1. Introduction 

Pulsars are obliquely rotating, magnetized neutron stars, with the obliqueness described 
by the angle, a, between the rotation axis and the magnetic dipole axis. The period, P, 
and the period derivative, P, are known for about 2000 radio pulsars. On a P-P diagram 
(actually a logP-logP diagram), these fall into three classes: normal pulsars, recycled pul- 
sars with smaller P, P, and magnetars with larger P, P. The conventional theory of pulsar 
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electro dynamics is built around two incompatible models, which we refer to as the vacuum- 
dipole model and the corotating-magnetosphere model. In the vacuum-dipole model the 
magnetosphere is ignored. Energy and angular momentum are carried away from the pulsar 
by magnetic dipole radiation at the rotation frequency, co = 27r/P. By equating the power 
radiated to the rate of loss of rotational energy, the model is used to derive the age, -P/2P, 
and surface magnetic field, Bsina oc {Ppy^"^, conventionally plotted as straight lines on a 
P P diagram. In the corotating-magnetosphere model, the magnetosphere is assumed to 
be corotating with the star, requiring the presence of the corotation electric field, whose 
divergence implies the Goldreich- Julian charge density, pcj- An additional simplifying as- 
sumption is made to reduce the electrodynamics essentially to electrostatics: an explicit 
assumption that achieves this is that the rotation and magnetic axes are aligned, sin a = 0. 
Alternatively this simplification is achieved through the assuniptiori that the magnetosphere 
is time- independent in a corotating frame (jScharlemann et allll978l ). The energy and angu- 
lar momentum are assumed to be carried away by a pulsar wind that results from plasma 
escaping along the open field lines in the polar-cap region, defined by those dipolar field lines 
that extend beyond the light cylinder radius, r\c = Pc/27i. This quasi-electrostatic, quasi- 
stationary model ties many of the details to the stellar surface, notably space-charge-limited 
flow, g aps, primary pa r ticles, pair formation front and the carousel model for drifting sub- 



pulses (jSturrock 



197ll : iRuderman fc Sutherland! 1 19751 : iFilippenko fc RadhakrishnanI Il982 



Cheng et al.lll986l ). There is no formal justification for the dichotomy between these two 
models, where one includes time-dependent fields only when the plasma is ignored, and the 
other neglects the time-dependent fields when the plasma is included. There is a nonzero 
displacement cur rent in an obliquely corot ating magnetosphere, but this is neglected through 
the assumption (JScharlemann et allll978l ) that the magnetosphere is tim e-independen t in a 
corotating frame. This and other criticisms of the existing paradigm (JMichell |200J) sug- 
gest that the formal basis of the theory needs to reconsidered. Indeed when intrinsically 
time-dependent electric fields are allowed, through inclusion of the displacement current in 
Maxwell's equation, the system is found to be violently unstable to the development of large 



amplitude oscillations tha.t lead to periodic bur sts of pair creation ( iLevinson et al.l 12005 



Beloborogov fc ThompsonI 120071 : iTimokhiru |2010| ). Moreover, the existing paradigm is not 
proving effective as an interpretative or predictive tool. This is especially the case for some 
recently identified, intrinsically time-dependent pulsar phenomena, not ably a link between 



nulling, mode switching, subpulse drift ing and abrupt changes in P (JKramer et al.l 12006 
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Weltevrede et al.l 120071 : iLyne et al.l 120101 ) . These phenomena seem to require a purely mag- 
netospheric interpretation, and the quasi-electrostatic, quasi-stationary theory tied to the 
stellar surface has had at best limited success as a basis for their interpretation. 



In this paper we consider the implications of the inductive electric field in an obliquely 
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rotating magnetosphere. The "inductive" electric field, Eind, which has curl Eind 7^ and div 
Eind = 0, results from the time- varying magnetic field, which is an essential ingredient in the 
vacuum-dipole model. Eind is ignored in the corotating-magnetosphere model. This neglect 
is implicit in the assumption that the magnetosphere is corotating, which implies that the 
only electric field in the magnetosphere is the corotation electric field, Ecor- However, even 
if corotating plasma is present, Ejnd is still generated by the changing magnetic field. The 
neglect of Eind is justified only if it is screened by the magnetospheric plasma. We point out 
that Eind cannot be screened by charges. In principle, the associated displacement current, 
£o<9Eind/c?t, can be screened by a current, and this might arguably screen the inductive 
electric field itself. We calculate the screening current density, and discuss whether or not 
current screening occurs in a pulsar magnetosphere. We conclude that complete screening of 
Eind cannot occur. An unavoidable consequence is that the electric field in the magnetosphere 
is not equal to E^or, so that the magnetosphere cannot be corotating. Inclusion of Eind in 
the theory changes the way we need to think about pulsar electrodynamics. 

In ^we calculate the fields in the vacuum-dipole model, first approximating it by a point 
dipole, and then generalizing to include the surface charge on a star (the Deutsch model). 
We write down the fields in the corotating-magnetosphere model for an oblique rotator. In 
[|3]we introduce the concept of current screening, and discuss the parallel and perpendicular 
components of the screening current separately. Perpendicular current screening is argued 
to be ineffective in a pulsar magnetosphere, and we discuss the resulting departure from 
corotation in §11 Possible implications for the interpretation of drifting subpulses and high- 
energy emission are discussed in §3 



2. Electric fields 

Exact expressions for the pulsar electric field may be written down for three mo dels: the 



vacuu m-dipole model for a point dipole, the Deutsch model for a centered dipole ((Deutsch 



19551 ). and an obliquely corotating magnetosphere. 



2.1. Vacuum-dipole model 

The fields from a time- dependent magnetic dipole, m(t), depend on the retarded time 
tret = t — r/c. The vector potential is 



A(t,x) = ^cur/^^^^^*)^-^° 
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where x is the position vector with respect to the center of the pulsar, and where a dot 
denotes a time derivative. The electric and magnetic fields are determined by 



E(t,x) 
The electric field is given by 



5A(t,x) 



B(t,x) =curlA(t,x). 



^('-)-£ 
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(2) 



(3) 



The magnetic field is given by 

3x X ■ m — r^m 3x x ■ lii — r^m x x (x x m) 



B(t,x) 






r^c 



ry"j z^-" 



(4) 



We refer to the term in equation (j4]) proportional to m as the dipole field, and the terms 
in equations ([3]) and dl]) proportional to rii and m as the inductive and radiative terms, 
respectively. In the following discussion, the radiative terms are ignored, except where 
stated otherwise; for most purposes they can be combined with the inductive terms. 



The magnetic field depends on time, with 
9B(t,x) 



dt 



-curlE(t,x) = ^ 

^ ' Air 



3x X ■ lii — r^m 3x x ■ m — r^rii x x fx x rii) 
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and it has a nonzero curl: 



,„, , l9E(t,x) 
curlB(t,x) = — ^^^ 



dt 
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2.2. Inductive electric field for a rotating dipole 

For a dipole rotating with angular velocity a;, one has 

rii = u; X [u; X (u; X m)]. 



m = a; X m. 



m = u; X u^ X m 



(7) 



Of particular interest in the following are the parallel and perpendicular components of 
the inductive field, where by 'parallel' we mean along the dipolar field lines. The parallel 
component is 



^indll (t, X) 



Eind(t,x)-B(t,x) 
|B(t,x)| 



/io X ■ mu) ■ m — X ■ w |m| 

'47rr2 [3(x-m)2 + r2|m|2]i/2 



(8) 
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where 'ind' denotes the inductive field. The time derivative of the parallel component of the 

inductive field is 

^p (f^\ = /^o X ■ (g; X m) u; ■ m 

^^i^iM|jir,xj 4^^2[3(x.ni)2 + r2|m|2]V2- ^^^ 

The perpendicular component of the inductive field implies a drift velocity given by 

_ Eind(t, x) X B(t, x) _ X ■ m[2r^Lj x m - 3xx ■ (a; x m)] 
^~ |B(t,x)|2 ~ 3(x-m)2 + r2|m|2 ' ^ ' 

where only the dipole field is retained in B(t, x). 



2.3. Deutsch model for a centered dipole 

The Deutsch model for a rotating magnetic star has the magnetic field inside a per- 
fectly conducting star described in terms of two functions. For a centered dipole, Deutsch's 
functions are Ri{r) = /iom/27rr^ and R2{r) = fiQm/ATrr^. The important change from the 
point-dipole model is the inclusion of the corotation electric field inside the perfectly con- 
ducting star. At the stellar surface (r = i?*) with a surrounding vacuum, the tangential 
component of the electric field must be continuous. This implies a surface charge density on 
the star, and a potential electric field outside the star. In an aligned model, this potential 
field is the only electric field present outside the star. In an oblique rotator, the Deutsch 
model has two contributions to the electric field in the vacuum outside the star: the poten- 
tial field (which depends on time) and the inductive field (unchanged from its point-dipole 
value). 

A potential field may be expanded in multipoles, and for a dipolar magnetic field, the 
potential electric field is quadrupolar. Its explicit form is 

o d2 

Equad(i^)X) = ^ [Sx w -xm-x — r^(xa;-m + ma;-x + a;m-x)]. (11) 

The particular form of Equad for an aligned rotator is written down in equation (ITSll below. 
The potential field plays a central role in the conventional model of pulsar electrodynamics. 
It is assumed to be screened by charges drawn from the stellar surface, with this screening 
breaking down in gaps, where the component of Equad along the magnetic field lines acceler- 
ates charges to high energy. The potential field plays essentially no role in the discussion in 
the present paper, and it is included here to emphasize that it is unrelated to the inductive 
electric field. 
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2.4. Corotating magnetosphere 

The corotation electric field is determined by the condition that there be no electric 
field in the local rest frame of the plasma, implying 

Eeor(t, X) = -(a; X X) X B(t, x). (12) 

We emphasize that the assumption of corotation requires that the electric field in the mag- 
netosphere be Ecor, and that if the electric field is not equal to E^or, the magnetosphere is 
not corotating. The following remarks relate to implications of the assumption of corotation. 



Any rotating vector field, V(t,x), satisfies an equation of motion (JMelrosd 119671 ) 

dVd x) 

)-^—^ = a; X V(t,x) - (w X x) •gradV(t,x) 

at 

= curl[(a;x x) X V(t,x)] - (a;x x)divV(t,x), (13) 

and Ecor(^, x) satisfies this equation. All electromagnetic fields in a corotating magnetosphere 
must satisfy equation (fT3l) . 

In an oblique rotator, the divergence, curl and time derivatives of Ecor are all nonzero. 
The divergence of Ecor(t,yi) determines the corotation charge density: 

divEcor(t,x) =p(t,x)/eo = -2a; • B(t,x) + (w x x) ■curlB(t,x). (14) 



For a dipolar field, equation ( 1T^ reduces to the Goldreich- Julian charge density (jGoldreich fc Julian 



19691 ). The divergence and time derivative (its displacement current) of Ecor are zero for an 



aligned rotator, a nd are specifically negle cted by the assumption of time-independent in a 



corotating frame (jScharlemann et allll978[ ). 



2.5. Electric fields in gaps 

To emphasize one role of the inductive electric field we discuss conventional models for 
gaps in an aligned rotator, and comment on how the parallel component of the inductive 
field requires a change in the interpretation. 

For an aligned rotator, the potential field can be written in the form 

Equad = -V^quad, "^quad = ~^^^ P2icOS 6) , (15) 

where $quad is the potential associated with the surface charge on the star, and where P2 is a 
Legendre polynomial. This vacuum field has a nonzero parallel component E\\ = Equad' B/|B| 
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along the magnetic field. In a corotating magnetosphere, the electric field given by equation 
f fTSjl is perfectly screened, and the only electric field is the corotation field. The concept of 
a gap is associated with charge starvation, which refers to the situation when there is an 
insufficient number of charges to provide the Goldreich- Julian charge density. A counterpart 
of equation (TT5|) then develops in the gap; the (potential) electric field in the gap can be 
attributed to surface charge densities on the lower and upper sides of the gap. Corotation is 
affected by a gap: the angular velocity of rotation changes across the gap. 



In an inner gap model (JRuderman fc Sutherland! Il975l ). charges from the stellar surface 
are assumed to provide the screening charge density immediately above the surface. In the 
absence of another source of charge, there is a deficiency in the number of charges needed to 
maintain pQj, and this deficiency increases with height. This leads to E\\ ^ developing in 
the inner gap; this E\^ accelerates primary particles to high energy in the gap, such that they 
emit 7 photons that produce secondary pairs. Charge separation between the secondary pairs 
provide s the additional charge density needed to screen E\\ at greater heights. In the original 
model (JRuderman fc Sutherlandlll975l ). the polar cap region above the gap sub-rotates, i.e., 
slower than the star to an observer, with E'y proportional to the difference in angular speeds 
across the gap, and determined by equation (TT^ with uj replaced by this difference and R 
interpreted as the radial distance to the gap. A criticism of the concept of a quasi-stationary 
gap is that when time-dependence is included through the displacement current, the gap is 
violently unstable to the development of large-amplitude electric oscillations, as discussed 
further below. 

Screening can also break down in other regions of the polar cap region, with the outer 
gap being a notable example. The Goldreich- Julian charge density is zero on surfaces where 
u; ■ B = 0, which corresponds to cos^ Q = 1/3, and has opposite signs on either side of this 
surface along any field line that passes through this surface. An additional source of charge 
is needed to allow this change in sign to occur. In the absence of an additional source of 
charge, screening is incomplete and the potential field redevelops , leading to E \ \ ^ 0. As 
in the inner gap, charges are accelerated and emit gamma rays (jCheng et al.lll986l ). The 
locatior i and details of an o uter gap model are affected by the current flow through the gap 
region ( iHinotori et al.ll2003l ). A simple model for the parallel electric fleld, E\\, in the gap is 



dE» 



£o- 



ds 



p = PGJ-e[N+{s) + N^{s)\ 



(16) 



where s denotes distance along the fleld line. The term in square bracket is the charge 
depletion from the Goldreich- Julian value. This electric fleld accelerates charges, leading to 
pair production. An outer gap is a favored location for the emission of the observed pulsed 
gamma rays from some pulsars. 



The electric fields in gaps in a conventional model are potential fields: they are caused 
by charges and can be screened by charges. In an oblique rotator, the presence of Ein^ 
implies an unrelated contribution to -Ey, from the parallel component of Eind- This inductive 
E^\ has a different functional form from the potential field; one has -Equad oc w/r^ and i^ind oc 
(cj/r^) sin a. Although relatively unimportant in the inner magnetosphere, E'ind is likely to 
be the dominant field in an outer gap region. 



3. Current screening 

In this section we discuss the concept of current screening in an oblique rotator. The 
inductive electric field is separated into perpendicular and parallel component, and only the 
parallel component can be screened by charges. In principle, the perpendicular component 
can be screened by a current. This current is identified, and whether or not current screening 
actually occurs is then discussed. 



3.1. Charge screening of the parallel inductive field 

The parallel component of the inductive field can be screened by charges. With s the 
distance along the field line, one separates the inductive field into parallel and perpendicular 
components, and writes 

div EiM = div^Eind^ + —^ = 0. (17) 

One can then identify a charge density p = — eo<9i?ind||/<9s. If this charge density is present 
in the plasma, it produces an E\\ that is equal and opposite to -Eindy, effectively screening 
i?ind||- Provided there is an adequate supply of charge, this screening should occur. Charge 
starvation can have a similar effect to that in conventional gap models: if there are insuffi- 
cient charges to screen i?ind|| completely, its presence leads to acceleration of particles, and 
associated pair production, until there are sufficient charges to restore the screening. 

This argument suggest that the role of Eind might be closely analogous to that of the 
potential field in an aligned model. However, there is a major difference when the displace- 
ment current is taken into account. Analytic and numerical solutions of the ID version of 
equation flT7|) show that it i s violently unstable to the build up of large- amplitude oscillations 



in th e parallel electric field (JLevinson et al.ll2005l : iBeloborogov fc Thompsorul2007l : iTimokhin 



20101). Screening of Eind|| can occur only in a time-averaged sense, where the average is over 



these oscillations. Pair creation in gaps in the conventional model is replaced by pair creation 
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at phases of the large amphtude oscillations. The concept of a spatially localized gap is no 
longer relevant. 



3.2. Displacement current 

The Maxwell equation 

1 9E(t,x) 



curlB(t,x) = yUoJ(t, x) + 



.2 



dt 



must always be satisfied. Current screening corresponds to the displacement current, which 
is the final term in ( !T8|) . being partially or completely replaced by the current J carried by 
the charged particles in the plasma. Before discussing screening, it is relevant to explain how 
equation ( 1T8|) is satisfied in the vacuum-dipole and corotating-magnetosphere models. 



For the inductive fields in vacuo, equation (ITSl) is satisfied with J = 0. The inclusion of 

a magnetosphere can change all three terms in equation flTSj) . We assume that the change 

to curl B is less important than the changes to the two terms on the right hand side. Ideal 

screening requires that the final term in equation ( ITSl) be zero. This determines the screening 

current density as 

Jsc(t, x) = -curl B(t, x) = , /Eind(t,x) ^ 
fiQ dt 

where Eind(i, x) is the inductive electric field that would be present in the absence of screen- 
ing. In identifying the screening current by equation (fTOll . we effectively require that the 
displacement current in vacuo be replaced by an identical J carried by charges. If the plasma 
cannot supply this current, then the displacement current is effectively unchanged from its 
value in vacuo, and Eind(i, x) must have essentially the same value as in vacuo. We ar- 
gue below that in a pulsar magnetosphere, the parallel component of equation flT9l) may be 
satisfied, but the perpendicular component cannot be satisfied. 

For the corotation field, the time-derivative of equation ( IT2|) can be evaluated using 
equation (fT3l) . giving 

"^ ' ' = curl [(a; X x) X Eeo,-(t, x)] - (a; X x) div E,o,(t, x). (20) 

One can rewrite equation ( l2Ql) in the form of equation (1181) . and re-interpret it. The interpre- 
tation is as a relation between the displacement current associated with Ecor(^,x), the curl 
of the corotation- induced magnetic field, Bcor(^, x) = (w x x) x Ecor(i, x)/c^, and the current 
density, pcj'^ ^ x, due to the corotating charge density. The corotation-induced magnetic 
field is smaller than B(t,x) by a factor of order r'^/rf^, and can be neglected in the inner 
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magnetosphere. Thus equation f lTSj) is satisfied by the corotation fields alone. This justifies 
the neglect of the corotation field in identifying the screening current given by equation flT9|l . 



3.3. Ideal screening in a corotating model 

The assumption that the magnetosphere is corotating has the implication that the in- 
ductive field is absent and hence must be screened. With the inductive field given by equation 
([3]), the required current density is 

X X [a; X (w X m)] , , 

The parallel component of the screening current follows from equation (jH]): 

J„rtx) = i x-(a;xm)u;-m 

We emphasize that the current density given by equation ( |2T]) is required by the hy- 
pothesis that the magnetosphere of an oblique rotator is rigidly corotating. In order for this 
to be the case, the plasma must supply the current Jgc, given by equation (IT9l) . 



3.4. Parallel current screening 

The parallel component of Eind can be screened by charges, and the parallel component 
of (9Eind/c^i can be screened by the current given by equation (l22l) . Such screening occurs 
provided that there are sufficient charges available to provide the charge and current densities. 
The required number of charges can be estimated as follows. 

The inductive and corotational fields are proportional to w/r^, and they differ in magni- 
tude only by geometric factors. Hence, the charge density required to screen the -EindU differs 
from the Goldreich- Julian value, pcj, only by a similar geometric factor. Screening of -Eind|| 
occurs provided tha t secondary pair creation results in a multiplicity, M, greater than unity 



( iBeskin et al.l 119931 ). Specifically, if the number densities of electrons and positrons are n±, 
and pgj = e(n+ — nS) is the number density required for corotation, then the requirement 
is that the multiplicity, M = (n+ + n_)/pQj, be greater than unity. 

The requirement on the parallel current density for screening of dEij^d\\/dt can be esti- 
mated by noting that the maximum current density is when the electrons and positrons are 
fiowing in opposite directions at relativistic speeds. This maximum is e(n_(_ + ra_)c = MpQjc. 
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Apart from factors of order unity, one has Jsc\\/MpQjc ~ r/Mric, which is much less than 
unity for r <^ r\c = c/u. It follows that pair creation (in large-amplitude oscillations) pro- 
vides an adequate source of charges to ensure that parallel current screening occurs in the 
inner magnetosphere. However, parallel current screening must break down in the outer 
magnetosphere. 



3.5. Perpendicular current screening 

Perpendicular current screening involves different physics from parallel current screen- 
ing. We discuss this from two complementary viewpoints. 

The equivalent dielectric tensor of any plasma at very low frequencies may be approxi- 
mated by a perpendicular component 1 + c^ /v\ and a parallel component 1 — bJijiJ^ ^ where 
the unit terms corresponds to vacuum. Here the Alfven speed is va = B/(jj,Qr]Y^'^, where t] 
is the mass density. In a conventional plasma, one has va <^ c, but in a pulsar plasma, one 
has Va ^ c. Hence, the perpendicular response of a pulsar plasma is effectively the same 
as if the plasma were absent. The perpendicular response links the current density, J_l, to 
the displacement current, 6odE±/dt. It follows that the plasma can supply only a fraction 
(?/v\ <^ 1 of the current required to screen the displacement current. In contrast, interpret- 
ing —oj"^ as a second time derivative, the parallel response correspond s to dJ \\ /dt = cu^gn-E' ii i 



which leads to large amplitude electric oscillations in a pulsar plasma ( jLevinson et al.ll2005l ). 



The perpendicular plasma response (at very low frequencies) can be understood in terms 
of the so-called polarization drift. Due to 9Ex/(9t, a particle with charge q and mass m drifts 
across the magnetic field lines at a velocity {'m/qB'^)dFjj_/dt. For an inductive electric field 
due to the magnetic field varying at frequency u = 2tt/P, the polarization drift is smaller 
than the drift caused by the inductive electric field by a factor of u/Qe (see equation (10)), 
where ^e = eB/m is the cyclotron frequency. Summing over all charges, this leads to a 
current density J_l = {c^ /v'j^eQd¥i^/dt, which reproduces the result imphed by the plasma 
response tensor, providing a physical interpretation of this response. 

In summary, the parallel component of the displacement current can be screened by a 
plasma current, which, however, is unstable to large-amplitude oscillations. The perpendic- 
ular component of the displacement current is essentially unchanged from its value in vacuo. 
It follows that the neglect of the perpendicular component of the inductive electric field in 
models for pulsar electrodynamics is not justifiable. Implications of including the inductive 
electric field are discussed in the remainder of this paper. 
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4. Departure of corotation 



The neglect of the inductive electric field in an obliquely rotating pulsar is not justified, 
and its inclusion implies that the magnetosphere cannot be in rigid rotation. 



4.1. Inductively induced drift velocity 

The perpendicular component of the inductive electric field implies an electric drift Av, 
determined by equation flTO|) . It is convenient to introduce spherical polar coordinates, r, 6*, cj) 
defined by the rotation axis. The spherical polar components of the dipolar field are 

(2 [cos a cos 6 + sin a sin 9 cos(0 — ut)] \ 
cosQ;sin6' — sinacos6'cos(0 — wt) , (23) 

sin a sin(0 — ujt) J 

where the initial conditions are chosen such that the magnetic axis is in the plane = at 
t = 0. The angle 6*^, defined by cos6'm = cos a cos 6^ + sin a sin 6* cos (0 — wt), corresponds 
to the magnetic colatitude, which varies periodically as the star rotates. At the phases 
(f) — Lot = nvr, n = 0, ±1, ■ ■ ■ , the field line is in an azimuthal plane, but at other phases the 
field has a nonzero azimuthal component {B^ ^ 0). 




The drift velocity given by equation ( jTOj) has spherical polar components 

„ , sin 6^ smid) — ut) 
ujr sma cos Um I „ \ \ , , •. 

2cos9sm{(P-ujt) . (24) 




I + 3C0S'^ 6m \ n ^J. ^\ 

2 cos(0 — Lot) 

The drift velocity is perpendicular to the magnetic field, implying AvrBr + AveBg + Av^B^ = 
0. The drift velocity has its extrema at cos(0 — Lot) = ±1, with values 

2a;r sina cos(a; =F 6') ,„^, 

The maximum drift speed occurs for cos(a =F 6*) = 1, and has value At>max = ^Lorsina. 



4.2. Small obliquity approximation 

In the case sin a -C 1, the effects of obliquity can be treated as a perturbation on the 
aligned model, as discussed in the Appendix. To zeroth order in sin a, the magnetosphere 
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is assumed to be rigidly corotating, with the charge density equal to pcj- The inductive 
electric field is of first order in sin a, 

Homu sin a / /^ .x 1 /or\ 

\ cos6' sin(0 — cut) ) 

and is unchanged to lowest order in the iteration. 

On expanding the drift velocity given by equation f|24p in sin a, the first order terms 
average to zero over a rotation period, and we retain a term of second order in sin a to show 
that this is not the case in general. The resulting expression for the drift velocity is 





bjr sin a cos 6* 
1 + 3 cos2 Q 



1 — 3 COS^ Q . . n / I \ 

1 H ;r— sm a sm cosf — ut) 

1 + 3 cos^ 6 



sin 9 sin(0 — ut) 

2cos9 sin(0 — ut) 

2 cos(0 — ut) 

(27) 
with the expression inside the square brackets replaced by unity to lowest order. For sin a <^ 
1, we have cos6'm = cos(q;=f6') ^ cos^ (see figure [1]). The average drift over a rotation period 
is in the direction and is of order sin^ a: 

/A V > . ^ > .9 cos6' (1 — 3cos^ 6^) _„, 

{A„,) = C-=.n9. C = Wa^^i:^^^. (28) 

This average drift vanishes in the equatorial plane, and has opposite signs in opposite (ro- 
tational) hemispheres. 

The drift given by equation ^F!\ is superimposed on the assumed rigid corotation. The 
drift does not correspond to a modification of the angular velocity. To see this consider a 
perturbation, 5uj, in the angular speed of the magnetosphere. This would give a perturbation 
5uj r sin 6 in azimuthal velocity; this has a maximum in the equatorial plane, has the same sign 
in both hemispheres and is independent of rotational phase. The velocity given by equation 
f l7r|) satisfies none of these conditions. The magnetosphere cannot be in rigid rotation at 
any angular velocity. 

Figures [21 [3] and H] show plots of Af,^, Avg and Avr, respectively, in units of ur as 
functions of ut with t = at = 0. Each curve represents values for 6 = 0°, 10°, 20° or 
30° for a = 30°. The correction in square brackets in equation (127]) is included, so that the 
temporal variations include a term varying sinusoidally as (p — ut with a correction term 
that includes a variation as 2(0 — ut). The components Av^ and Avr are odd functions of 
cos 6', implying that they have opposite signs in opposite hemispheres, whereas Avg is an 
even function of cos 9. 



-14- 



The pulse window corresponds to a small range of </> and 9, and because an observer 
can see the drift only within this window, only a correspondingly small range in figures [2l 
12] and m is relevant to observations. This range is model-dependent, and in most models 
it corresponds to a range Acf) about </> = 0° and a range A^, of order A^, about a line of 
sight, 9, close to a. The relevant regions in figures |2l [3] and H] are near wt = for the solid 
curves. The three curves for 9 = a = 30° are replotted in figure [51 As the variation with 9 
is relatively small, the observable (relative to corotation) drifts correspond to the range A0 
about wt = in figure O The components Af ^ and /S.Vr are zero at uot = 0, implying that 
the observable drift is in the direction. The magnitude of this drift may be estimated by 
setting 9 = a,(f) — ut = in equation f l271) . This implies an observable drift 



{^Vcf>)ohs 



2ujr sin a cos a 
1 + 3 cos^ a 



1 + 



1 — 3 cos^ a 
1 + 3 cos^ a 



sin^a 



(29) 



Note that the sign of (Af<^)obs depends on the sign of cos a, that is, on the sign of u;- B. 
This sign also determines the sign of the Goldreich- Julian charge density. The inductively 
induced drift implies super-rotation when the Goldreich- Juhan density corresponds to an 
excess of electrons, and to sub-rotation when the Goldreich- Julian density corresponds to an 
excess of positrons. 



5. Discussion 

The inductive electric field, Eind, (due to the obliquely rotating dipole) cannot be com- 
pletely screened in a pulsar magnetosphere. Its presence implies a plasma drift across the 
field lines, such that the motion of the magnetosphere is not rigid corotation with the star. 
We explore the suggestion that this drift might provide a natural explanation for subpulse 
drifting. We also comment on the possible role of the inductive field in an outer-gap model 
for pulsed gamma-ray emission. 



5.1. Subpulse drifting 



The inductively-induced drift relative to corotation suggests a natural explanation for 
subpulse drifting. Drifting subpulses correspond to a systematic motion of the plasma from 
which the radio emission escapes within the pulse window corresponding to a small range 
of phases during which emission is obse rved. Data on subpulse drifting in a large sample 
of pulsars (IWeltevrede et al.ll2006l . 120071 ) lead to several general properties. One is that the 
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subpulse drift rate is independent of frequency. This is explained naturally by this drifting 
model, which is essentially geometric and varies only slowly with height, and hence with 
frequency assuming frequency-to-radius mapping. Another property is that the number 
of pulsars with positive and negative drift rates are roughly equal. This has a natural 
explanation due to the drift rate depending on the sign of cos 9. For small impact parameter, 
implying 9 k. am the polar cap region, the drift rate has one sign for pulsars with a < 7r/2 
and the opposite sign for pulsars with a > tt/2. The sign of the drift rate depends on the 
sign of the Goldreich- Julian density above the polar caps. 

A simple prediction of the model is that the subpulse drift rate should be proportional 
to sin a. This follows from the fact that the drift rate is proportional to sin a for small sin a. 

A more subtle prediction is based on the sign of the drift rate depending on the sign 
of cos 6', which is interpreted as the sign of cos a for emission from the polar cap region. 
The sign of cos a determines the sign of pcj • Some pulsars have a pr eferred handedness fo r 



circular polarization, and this is plausibly related to the sign of pcj (JMelrose fc Ludl2004j ). 
The prediction is that the sign of the subpulse drift rate and the handedness of the net 
circular polarization should be correlated. 



5.2. Acceleration by the inductive electric field 

Pulsed high-energy emission from pulsars is attributed to acceleration of particles in 
the outer magnetos phere, for example, due to the breakdown of screening in an outer gap 



( iTakata et al.ll2010l ). In such models, the parallel electric field in the outer gap is assumed to 
be a potential field. The inductive electric field, which is not included in any existing model, 
can play the same role as the assumed potential field, suggesting a natural alternative to 
outer-gap models. Although the inductive field calculated in section |2] becomes invalid as 
the light cylinder is approached, it can be used to provide a rough estimate of the actual 
electric field. This estimate is E\\ ^ uB^Rl/r"^ ^ Bicc{ric/r)^, where Bi^ = B^^{R^,/ric)^ is the 
magnetic field at the light cylinder. The main difference between such an inductive model 
and the conventional outer-gap models is that the parallel electric field results from the 
time-varying (obliquely rotating) magnetic field, and is not due to a local charge separation. 
The parallel inductive electric field appears when there are too few charges to provide the 
charge density needed to screen it. 

There is a close analogy between screening of the inductive field in the magnetosphere 
(r < ric) and screening of the radiative field in the wind zone (r ^ r\c). Early in the 
development of pulsar theory it was recognized that the radiative component of the vacuum 
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field, given by equations (|3]) and (jl]), would accelerate test particles to extremely high energy 
( JGunn fc Ostrikerl Il969[ ). In the theory of the pulsar wind, developed subsequently, the 
electric field is specified by the wind equations, and it is implicit that the radiative field is 
perfectly screened. This screening, and the generation of the electric field in the wind model 
require a current in the wind. If t here are too few charges (charge or current starvation) 
in the wind to sustain this current (jUsovlll994l ). the vacuum field is incompletely screened 
and accelerates charges. The breakdown of perfect screening leads to a transfer of energy 
from the Poynting fiux to a kinetic energy fiux, providing a natura l expl anation for the 
dominance of the latter in the outer wind zone (JMelatos fc Melrosd Il996l ). Similarly, in 
the inner magnetosphere, an unscreened inductive field accelerates a test charge along a 
magnetic field line to very high energy. Breakdown of parallel current screening is expected 
in the outer magnetosphere (r < ric). The parallel component of the inductive electric field 
then reappears, and accelerates charges. This needs to be examined in detail as an alternativ e 
to an outer gap for the production of pulsed 7 rays (JHinotori et al.ll2003l : iTakata et al.ll2010l ). 



In a conventional pulsar wind model ( jRees fc GunnI Il975l : lAronsI |2004| ) the magnetic 
field well beyond the light cylinder has a predominantly toroidal component. The rotation 
of the star and the outfiow of the wind cause the direc tion of the toroid al component to 



reverse periodicall y with radial distance ( iKirk et al.ll2007l ). Specific models (lBogovalovlll999 



Contopoulosll2005[ ) imply a temporally changing electric field in the wind. In such models the 
fiuid theory determines the elec tric field. Althou gh this electric field is varying periodically 
with time in an oblique rotator (JBogovalovlll999l ). implying a nonzero displacement current, 
this is not the displacement current associated with the rotating dipole. The latter field 
would be present in vacuo, and its implicitly assumed absence in wind models requires that 
it be effectively screened by currents. As with the inductive field in the inner magnetosphere, 
if current screening is ineffective, the radiative electric field must be present in the wind zone. 



6. Conclusions 



The neglect of the inductive electric field, Ejnd, in models for a pulsar magnetosphere 
cannot be justified. This field is generated by the changing magnetic field and the associated 
displacement current of an obliquely rotating magnetic dipole. It is impossible in principle 
to screen an inductive field by charges. We show that screening of the displacement current 
is possible in principle, and we calculate the required screening current density in an ideal- 
ized model. Screening of the components parallel to the magnetic field and perpendicular 
to the field involves different physics. The parallel component is unstable to the develop- 
ment of oscillations, and screening can occur only in an average sense, where the average 
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is over the oscillations. Perpendicular screening involves a polarization current, driven by 
the displacement current, and can never be complete. Nearly complete screening is possible 
for v\l(? <^ 1, which is not the case in a pulsar magnetosphere, where the perpendicular 
inductive field is essentially unscreened, and unchanged from its vacuum value. 

The presence of Eind has important consequences for understanding pulsar magneto- 
spheres. In particular, the magnetospheric plasma cannot be corotating with the star. The 
assumption of corotation is central to conventional corotating-magnetosphere model for pul- 
sars, and much of the physical interpretation of pulsar phenomena is based on this assump- 
tion. The inductive electric field is absent only in the artificial case of an aligned rotator 
(sin a = 0). In the case of small obhquity, sin a ^ 1, an aligned rotator can be regarded 
as a zeroth order model, with the inductive field and its implications treated to first order 
in sin a. Using such a model we calculate the inductively induced drift rate, which is the 
difference between the actual velocity of the local magnetospheric plasma and the corota- 
tion velocity. We suggest that this drift offers a natural explanation of drifting subpulses. 
The drifting of subpulses is a direct result of the perpendicular component of the inductive 
electric field. This suggests that subpulse drift should be interpreted as a signature of the 
inductive electric field and used to infer properties of this field. 

The parallel component of the inductive field can be screened by charges, and the 
breakdown of this screening (due to charge starvation) in the outer magnetosphere offers 
an alternative to conventional (outer gap or slot gap) models for pulsed gamma ray emis- 
sion. Conventional models invoke E\\ due to separation of charges. The conceptual change 
suggested here is that the (parallel) inductive field would be present in the absence of the 
screening, and it must appear whenever charge or current starvation limits the ability of the 
plasma to screen it. 

A surprising implication of including the inductive electric field is that the magneto- 
sphere cannot be corotating with the star. The conventional argument is that MHD requires 
that the electric field be zero in the rest frame of the plasma, and hence is equal to the Ecor 
in the pulsar frame. Inclusion of the inductive electric field implies an inductively driven 
drift, and this drift motion is implicitly and incorrectly assumed to be zero in a corotation 
model. 

Finally, we endorse a remark that the neglect of the displac ement current in ast rophysi- 



cal and space plasmas has led to conceptual misunderstandings (jSang fc Lysakll2006l ). In the 
context of pulsars, the conventional quasi-electrostatic corotating-magnetosphere model ex- 
cludes the displacement current, but when the displacement current is included, the system 
i s found to be violently unstable to the developnaent of large amplitud e electric oscillations 



(JLevinson et al.l l2005l : iBeloborogov &: Thompson! 120071 : iTimokhin 



2010|) 
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A. Self-consistent fields 

Both the vacuum-dipole and the corotating-magnetosphere models are based on assump- 
tions that are not satisfied. The vacuum-dipole model neglects the plasma, but the plasma 
screens the parallel components of the inductive electric field and the displacement current, 
invalidating their derivation using the model. The corotating-magnetosphere model is based 
on the hypothesis that the magnetosphere is corotating, but the perpendicular component 
of the inductive electric field is essentially unscreened which implies an electric drift incom- 
patible with corotation. This raises the question as to how, in principle, one can formulate a 
self-consistent model for the electrodynamics of an obliquely rotating pulsar magnetosphere. 
In this appendix, some remarks are made on how this might be achieved. 



This preprint was prepared with the AAS IATJt;X macros v5.2. 
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The determination of the fields and induced charge and current densities requires solving 
a self-consistency problem. The potentials due to given charge and current densities satisfy 

A(t,x) ; 4,ri |x-x'| V AoJ(«',x') )' ^ ' 

with t' = t — |x — x'l/c. Self-consistency results from requiring that the current density be 
minus the component of the displacement current along the magnetic field lines, and that 
the charge density is determined by the divergence of the perpendicular component of the 
electric field. 

The solution of equation flAip requires boundary conditions at the star. As in the 
Deutsch model, the corotation field inside the star determines the electric field at the surface, 
and the self-consistent field must join onto the interior solution continuously across the 
surface. This implies that the plasma immediately above the surface is nearly corotating, 
and deviations from corotation increase with height above the surface. 

To describe the effects of partial current screening in detail one needs to solve equa- 
tion (lAip . The screening cancels only the parallel component of the displacement current, 
Eod'Eiad/dt, and we are interested in the resulting modified inductive field. 

One can solve equation flAip iteratively for sufficiently small sin a. Identifying the 
zeroth order solution as the unscreened field due to the rotating dipole in vacuo, the current 
density is identified as the parallel component of Jsc, given by equation fl2Tl) . and the current 
in equation flAip is identified as this current directed along the dipolar field lines. The 
charge density is determined by integrating dp/dt = dJsc\\/ds, where s denotes distance 
along the dipolar field line, with respect to time. With these terms on the right hand side of 
equation (]A1|) . the solution is used to calculate the first order correction to the electric field 
and the associated current and charge densities. The second order corrections are found by 
repeating the calculation with the first order charge and current densities. To lowest order, 
the inductively induced drift modifies the rigid corotation that applies to zeroth order. 

Hence the inductive electric field is given by, 

EUt, x) = E^) (t, x) + E;;) (t, x) + ... + E;2(t, x) + E[^^'\t, x) + ..., (A2) 

where in the Lorentz gauge, E^^^^ (t, x) is a function of Ej\j^(t, x), which can be determined 
from 0(^)(t, x) and A('')(t, x) with, 

A(^+i)(t,x) y* " 47r y "^ "" |x - x'l V f^oJ^'Kt',^') J ' ^^'^> 
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where E-\j^(t, x) represents the unscreened field due to the rotating dipole in vacuo. 

The model used in section H] corresponds to the first order term in the iteration in 
sin a <^ 1. The iteration procedure must converge rapidly if sin a is sufficiently small, but 
there is no simple way of determining the range of validity of the approximation to first order 
in sin a. 
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Fig. 1. — The polar angles used in this paper are shown in the plane that contains the rotation 
and magnetic axes. Angle 9 is the polar angle relative to the rotation axis, and we denote 
the magnetic colatitude by 6m- This differs from a notation used in the pulsar literature 
( jLyne fc Manchester! Il988l ) where these two angles are denoted C, and 9, respectively. The 
angle (3 = 9 — 9m is referred to as the impact parameter. For sin a <^ 1, we have a ^ 9 ^ 9m- 
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Fig. 2. — The azimuthal drift velocity, /S.v^{t,0), as given by equation ( !27l) . plotted against 
time or rotational phase for a = 30°. The different curves represent different 6 values 
range from 0° — 30°. In addition to varying sinusoidally in magnitude, Av<^(t, 9) also varies 
asymmetrically with the positive maximum value different from the negative minimum value, 
which indicates that the drift velocity is faster in the second half of the rotation. 
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Fig. 3. — As for figure E] but for the polar component, Av0{t,6), as given by equation ( 1271) . 
The amphtude of a curve decreases as 6 increases as oppose to the other two components. 
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Fig. 4. — As for figure |2]but for the radial component, At>r(t, 6'), as given by equation ( l27jl 
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Fig. 5. — The three components Af^, Afg and Af^ are shown in the case when the hne of 
sight coincides with the magnetic axis at ^ = 30° 



